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ABSTRACT 

We use weak gravitational lensing to measure the masses of five galaxy clusters selected from the 
South Pole Telescope (SPT) survey, with the primary goal of comparing these with the SPT Sunyaev- 
Zel'dovich (SZ) and X-ray based mass estimates. The clusters span redshifts 0.28 < z < 0.43 and 
have masses M500 > 2 x 10 14 h~ 1 MQ, and three of the five clusters were discovered by the SPT 
survey. We observed the clusters in the g'r'i' passbands with the Megacam imager on the Magellan 
Clay 6.5m telescope. We measure a mean ratio of weak lensing (WL) aperture masses to inferred 
aperture masses from the SZ data, both within an aperture of i?soo,SZ derived from the SZ mass, of 
1.04 ± 0.18. We measure a mean ratio of spherical WL masses evaluated at i?5oo,sz to spherical SZ 
masses of 1.07 ± 0.18, and a mean ratio of spherical WL masses evaluated at i?5oo,WL to spherical 
SZ masses of 1.10 ± 0.24. We explore potential sources of systematic error in the mass comparisons 
and conclude that all are subdominant to the statistical uncertainty, with dominant terms being 
cluster concentration uncertainty and iV-body simulation calibration bias. Expanding the sample of 
SPT clusters with WL observations has the potential to significantly improve the SPT cluster mass 
calibration and the resulting cosmological constraints from the SPT cluster survey. These are the first 
WL detections using Megacam on the Magellan Clay telescope. 
Subject headings: cosmology: observations - galaxies: clusters: individual 
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1. INTRODUCTION 



The abundance of galaxy clusters as a function of mass 
and redshift is sensitive to dark energy and other cosmo- 
logical parameters through the growth function of large- 
scale struc ture (LSS) and the cosmological volume ele 
ment (e.g. , | Wang & Steinhardt|19 98 ; Ha iman et al.|2001 



Holder et al.||2UUl| |Weller fe Ba ttyc 200 3)) . As empha- 



sized by the Dark Energy 'i'ask Force (Albrecht et al 



2006 1 , the abundance of clusters provides constraints on 
dark energy that are complementary to those of distance- 
redshift relations, such as standard candles and rulers 
including Type la supernovae (SNe) and baryon acous- 
tic oscillations (BAO). Recent results using this method 
have shown that cluster surveys can significantly improve 



particularly the dark energy equation of state (|Vikhlinin 


et al. 


2009b; Mantz et al. 2010b 


JRozo et al. 2010; Benson 


et al. 

1 1 M_ 


2011||Reichardt et al.l|201 





to search for high-redshift, massive clusters, whi ch are 



particularly usefu l for constraining cosmology (e.g., Carl 



strom et al. 



scattering o 



2002). The SZ effect is the inverse Compton 
cosmic microwave background (CMB) pho- 
tons by the hot electrons in the intracluster medium. SZ 
observables are nearly redshift independent, and more- 
over, arc expected from simulations and observations to 
trace total cluster mass with low intrin sic scatter (e.g., 
Kravtsov et al.|2006)jBenson et al.|2011 ). To extract con- 
straints on cosmological parameters, the cluster redshifts 
must be measured with optical and infrared follow-up 
observations and the cluster masses must be estimated 
using accurately calibrated proxies. 



The South Pole Telescope (SPT; |Carlstrom et al.|20lT | 
is a millimeter-wavelength telescope that recently com- 
pleted a 2500 deg 2 SZ cluster survey. A catalog from 
the first 720 deg 2 of the survey has been released that 
includes 224 cluster candidates, 158 of which have con- 
firmed optical or infrared galaxy cluster counterparts 
with redshifts as high as z = 1.37, a median redshift 
of 0.55, and a median m ass of M^ ny « 2.3 x 10 14 /i _1 Af Q 
(jReichardt et al.||2012| hereafter |R12|). Using the 100 



cluster candidates at z > 0.3 above the SPT 95%-purity 
threshold, the SPT cluster data have been combined with 
CMB+BAO+iJo+SNe data to provide constraints on the 
dark energy equation of state of w = —1.010 ± 0.058, a 
factor of 1.3 improvement over the constraints without 
the cluster abundance data. However, this improvement 
was limited by the ~ 10% uncertainty in the SPT cluster 
mass calibration. 

The method to determine the masses of clusters in 
the SPT survey is described in detail by |Benson et al.| 
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(2011 hereafter Bll). In brief, the mass estimates are 
computed from the SPT SZ data, and are calibrated to 
scale with total mass using a combination of X-ray ob- 
servations and cluster simulations. For a subset of the 
SPT survey cluster sample, X-ray observations have been 
obtained to measure Yx, the product of the X-ray de- 
rived gas mass and core-excised temperature. This is 
used in combination with a Yx-mass relation that is cal- 
ibrated at low-redshift (z < 0.3) using X-ray d erived hy- 
drostatic ma ss estimates of relaxed clusters (Vikhlinin 
et al.||2009a). The calibration of the Yx~ mass relation 



is expected to be a ccurate to 10% based on simulations 
( Nagai et al.||2007 ), and has been empirically verified to 
have this level of accuracy usin g weak lensing (WL) ob- 
servations of the same clusters (iHoekstra 2007). 



WL here refers to the subtle tangential shearing of 
extended sources behind cluster halos, at projected dis- 
tances well outside the Einstein radius (for a review, see 
Bartelmann & Schneider 2001 ) . Gravitational lensing 
is sensitive to total projected mass and has the benefit 
of being insensitive to the dynamical state of the lens: 
the observables are independent of whether the cluster 
gas or galaxies are in hydrostatic equilibrium. WL mass 
estimates have been used to test the accuracy of X-ray 
mass estimates (e .g., Hoekstra 2007 Mahdavi et al. 2008 
Zhang et al.|2008 1; however, these works have not e xplic- 
itly cross-che cked the Yx-mass calibration of |Vikhlinin 



et al. (2009a 



and have used clusters predominantly at 
low redshift \z < 0.3) that are outside the SPT survey 
region. 

As the first steps toward a WL-based calibration of the 
SPT cluster sample, we have observed five clusters from 
the SPT survey with the Magellan Clay-Megacam CCD 
imager. The five clusters span redshifts 0.2 8 < z < 0.43 
and were selected from the SPT catalog of |R12| from a 
subset of clusters with existing or scheduled observations 
by either the Chandra or XMM-Newton X-ray space- 
telescopes. From this SZ-plus-X-ray sample, we ran- 
domly selected five clusters that were observable during 
awarded telescope time and which were at 0.3 < z < 0.5. 
The X-ray analyses for these clusters are in progress. We 
have also obtained multi-object spectroscopy of cluster 
members for precise redshift measurements of four of the 
clusters. 

The primary goal of this paper is to provide the first 
direct com pariso n of the SZ and X-ray derived mass es- 
timates of IR12I to WL mass measurements. We em- 
ploy aperture masses, which are largely insensitive to the 
properties of cluster cores, and have been shown in pre- 
vious works to s cale well with other, low-scatter mas s 
observables (e.g., |Hoekstra||2007| [Mahdavi et al.||2008 ). 
We also use spherical masses by fitting the WL shear 
data to analytic profiles. 

Magnitudes are in the Sloan Digital Sky Survey 
(SDSS) AB system unless otherwise noted. We adopt 
a flat ACDM cos mology with Qm = 0.27 and Hq = 
70kms _1 Mpc" 1 rtKomatsu et al.|201l| . Masses M A are 
defined at radius -Ka , where the mean interior density is 
A times the critical density of the universe at the cluster's 
epoch, pcrit(z) = 3H 2 (z)/8ttG, and H(z) is the Hubble 
parameter. 

2. THE CLUSTER SAMPLE 
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In this section, we provide an overview of the larger 
sample of clusters from which we selected targets for WL 
follow-up. We briefly discuss SZ detection with SPT, the 
SZ cluster center determination, and SZ mass estimation. 
We then discuss the cluster spectroscopic redshift mea- 
surements of four of the five systems selected for WL ob- 
servations. Table [T] summarizes these basic cluster data. 

2.1. SZ Cluster Detection 



The clusters were selected from the IR12I SPT cluster 
survey catalog. That work contains details of the sam- 
ple, the SPT data from which it was extracted, and the 
cluster extraction process. In summary, ^720 deg 2 of 
sky were surveyed by the SPT in the 2008 and 2009 ob- 
serving seasons to a depth such that the median mass 
of a cluster detected is M500 ~ 2.3 x l0 li h~ 1 M Q . Clus- 
ter candidates are extracte d from the data usi ng a mul- 
tifrequency matched filter (Melin et al. 2006). Twelve 
different matched filters are used spanning a range of 
angular scales for the assumed cluster profile, and the 
cluster candidates are ranked by the maximum detection 
significance across all filter scales, defined as £. All can- 
didates with £ > 4.5 are included in the catalog, with 
£ also used as the primary observable to determine the 
cluster mass. SZ significance maps for the five clusters 
discussed in this work are shown in the Appendix. 

We use the SZ detection positions as the cluster cen- 
ters in the baseline WL analysis, although we explore 
the effect of using other positions as well. The statistical 
uncertainty in SZ-determined positions is a function of 
the cluster size, the SPT beam size (FWHM = 1'.6 and 
1'19 at 95GHz and 150GHz respectively), and the signif- 
icance of detection. In the limit that clusters are point 
sources in the SPT data, the rms positional uncertainty 
is ~ 0fwhm/£, where #fwhm is the beam width ( [Ivi-| 
son et aL | [2007). For re solved clusters, the uncertainty 
is ~ V^fwhm + (k6 c ) 2 /£,, where 6 n is the cluster core 



size, and k i s a fa ctor of order unity Story et al. (2011); 
Song et al. (2012). For clusters such as those described 
in this work (with £ ~ 8 and z ~ 0.4), this uncertainty 
is estimated to be 10"— 15". 

2.2. SZ Mass Estimates 



In this work, we use the mass estimates from R12 



The 

SPT mass calibration and method to estimate cluster 
masse s from th e SZ and X-ray data is described in de- 
tail in |B11| and |R12| In summary, a probability density 
function of each cluster's mass estimate was calculated at 
each point in a Markov Chain Monte Carlo (MCMC) that 
varied both the cluster scaling relations and cosmological 
parameters assuming a ACDM cosmology and using the 
CMB+BAO+SNe+SPT CL data, where SPTc L denotes 
the added cluster data set. In effect, this step is calculat- 
ing the posterior probability given the measurement un- 
certainties and the expected distribution of galaxy cluster 
masses for that specific cosmology and scaling relation. 
The resulting masses are SZ plus X-ray posterior mass 
estimates where applicable; even for clusters without X- 
ray data, the mass normalization from clusters with X- 
ray data affects the SZ scaling relation parameters ex- 
plored by the chain. The probability density functions 
for different points in the chain are combined to obtain a 
mass estimate that has been fully marginalized over all 



cosmological and scaling relation parameters. The final 
products are labeled M 50 o,sz here. 

The uncertainty on Msoo.sz is dominated by a ~ 15% 
intrinsic scatter in the SZ-mass scaling relation, a ~ 10% 
uncertainty due to the finite detection significance of the 
cluster in the SZ maps, and a ~ 10% systematic uncer- 
tainty associated with the normalization of the SZ-mass 
scaling relation. Together, these and other sources of un- 
certainty yield a ~ 20% uncertainty on M 50 o,sz for each 
cluster. 

2.3. Spectroscopic Redshifts 

We have targeted four of the clusters for multi-object 
spectroscopy observations, the details of which will ap- 
pear in J. Ruel et al. (in preparation). In summary, we 
used GISMO (Gladders Image Slicing Multi-slit Option) 
on IMAC S (Inamori Magellan Areal Camera an d Spec- 



trograph; Dressier et al. 
Magellan Baade 6.5 m te 



|2003| |Osip et al.||2008[) at the 
descope in 2010 September and 
October. These observations used the //4 camera, the 
zl 430 — 675 filter, and the 300 lines mm -1 grating. Con- 
ditions were photometric and the seeing varied between 
075 and 0'.'9. Reduction of the raw spectra was done with 
the COSMOS package 35 (Carnegie Observatories System 
for Multi-Object Spectroscopy); redshifts were measured 
by cross-correlation with the "fabtemp97" template in 
RVSAO (Radial Vel ocity Smithsonian Astr ophysical Ob- 
servatory package, Kurtz & Mink 1998[ ) and checked 
for agreement with visually-identified features. Outliers 
were rejected by iterative clipping at 3a. The redshift of 
each clu ster was calculate d using the robust biweight es- 
timator ( Beers et al. 1990[ ) and the confidence interval by 
bootstrap resampling. These results are given in Table 

ED 

The redshift of SPT-CL J2030-5638 is photometrically 
derived from the red sequence of cluster galaxies , the 
details of which are described by Song et al. (2012) and 

ma 

3. WEAK LENSING DATA 

In this section we describe the acquisition, reduction, 
and calibration of the images and photometry from which 
weak lensing masses are measured. 

3.1. Observations 

We imaged through <?', r', and i' passbands using 
Megacam on the Magellan Clay 6.5 m telescope at La s 



Campanas Observatory, Chile (McLeod et al. 1998 



Megacam was previously commissioned on the Multiple 
Mirror Telescope, where it was used to study wea k lens- 
ing by galaxy clus ters in the northern hemisphere ( Israel 
et al.||2010||2011[ ). Megacam consists of a 9 x 4 CCD 
array producing a 25' x 25' field- of- view. We operated 
read-out in 2 x 2 binning mode for an effective pixel scale 
of / .'16. The r'-band images are used for shape measure- 
ments (Section^, and the added g' and i' bands are used 
in concert with deep photometric-redshift catalogs from 
external surveys to prune and characterize the source 
pop ulati on using magnitude and color information (Sec- 
tion [42 1. 

Observations of SPT-CL J0516-5430 occurred at the 
end of the nights of 2010 October 6 and 7. A total of 
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TABLE 1 

The cluster sample 



Cluster Name 




N 




SZ R.A. 


SZ Decl. 










(deg J2000) 


(deg J2000) 


SPT-CL J0516-5430 a 


0.294(1) 


48 


9.42 


79.1480 


-54.5062 


SPT-CL J2022-6323 b 


0.383(1) 


37 


6.58 


305.5235 


-63.3973 


SPT-CL J2030-5638 b 


0.40(4) 




5.47 


307.7067 


-56.6352 


SPT-CL J2032-5627 c 


0.284(1) 


31 


8.14 


308.0800 


-56.4557 


SPT-CL J2135-5726 b 


0.427(1) 


33 


10.43 


323.9158 


-57.4415 



Note. — Basic data for the five clusters we have targeted for weak lensing 
analysis. 

Column 1: R12 cluster designation. 

Column 2: mean cluster redshift as measured from the ensemble of cluster 
members. The number in parentheses is the uncertainty in the last digit. All 
redshifts are spectr oscopic except for SPT -CL J2030-5638, which was photo- 
metrically derived [Son£etaLj2012]| . 

Column 3: number ol cluster galaxies for which we successfully measured spec- 
troscopic redshifts. 

Column 4: peak signal-to-noise ratio of the SZ detection. 
Column 5: right ascension of the SZ signal-to-noise ratio centroid. 
Column 6: declination of the SZ signal-to-noise ratio centroid. 

a Discovered by |Abell et ah] | |1989[ l where it was design ated as AS0520. The 
spect roscopic redsnut was measured to be z; = 0.2950 | |Leccardi &: Molendi| 
[2008] , 



Not known prior to |R12[ 
c A number of st r uctu res have been identified near this cluster by other authors. 
See Section p7T2||R12| and |Song et al.| pol^ for discussions. 



1760 s of exposure time was obtained in the r'-band us- 
ing a 2 x 2 square dither pattern of 6" on a side, plus the 
same 2x2 pattern executed about 25" north and 56" 
west, for a total of eight individual exposures of 220 s 
each. The two star guiders, situated on opposite sides 
of the Megacam focal plane, were simultaneously op- 
erational for every r'-band exposure, which resulted in 
uniform and stable point-spread function (PSF) FWHM 
patterns across the entire field, as monitored upon read- 
out during observation. Twelve hundred seconds of total 
exposure time was obtained in the <?'-band, and 3600 s 
in the i'-band, each with dithers that covered the chip 
gaps. Conditions were clear and stable with good seeing. 

The rest of the clusters were observed over 3 second- 
half nights on 2011 May 31-June 2. For each cluster, a 
total of 1800 s of exposure time was obtained in the r'- 
band using a three-point diagonal linear dither pattern 
that covered the chip gaps. In the g' band, we obtained 
1200 s of exposures with the same three-point dither pat- 
tern, and in i 1 we integrated for 2400 s using a five-point 
diagonal linear dither pattern. Conditions in this run 
were intermittently cloudy, resulting in approximately 
50% unsuitable time. Seeing was sub-arcsecond in the 
r'-band for these four clusters. 

For all five clusters, special care was taken to observe 
in the r'-band in stable, good seeing conditions under the 
clearest available skies. We executed the dither patterns 
in immediate succession and monitored the seeing. The 
result was < 15% variation in r'-band seeing in all dithers 
for any given cluster. Because of this and the uniform 
PSF pattern afforded by Clay and Megacam (Section EJ), 
we coadd images without homogenizing to a common 
PSF in any way. 

3.2. Image Reductions 

The images are reduced at the Smithsonian Astro- 
physical Observatory (SAO) Telescope Data Center us- 
ing the SAO Megacam reduction pipeline. The pipeline 



uses publicly available u tilities from IRAF, SEx tractor 
( Bertin fc Arnouts|l996 ) , and Swarp programs 36 ( Bertin 



et al. 



, as well as in-house routines. 



Basic CCD processing includes overscan correction, 
trimming, and bad pixel removal. Cosmic r ays are re- 
moved using L. A. Cosmic (van Dokkuml 2001 ). Flat-field 



images for each filter were generated from sets of twi- 
light sky observations, taken during clear dawn or dusk 
twilight, and then applied to the data. 

To correct for any scattered light remaining after flat- 
fielding, an illumination correction is also performed. 
The illumination correction image for a given filter is 
made by using 12 exposures from a dithered pattern de- 
signed to expose the same stars over the full span of the 
focal plane. Where these data are not photometric, i.e., 
there is too much scatter in the residuals of the dithered 
exposures, we instead use just one exposure of the SPSS 
Stripe 82 field and match to the catalog of |Ivezic et al.| 
p007| . 

A fringe correction is performed for the i' filter. We 
make fringe-frames by combining large numbers of i' 
science-field exposures. 

In the final step, two passes are made of the world co- 
ordinate system (WCS). The first pass fits star positions 
to a refe rence catalog, the Tw o-Micron All Sky Survey 
(2MASS, [Skrutskie et al.|2006[ ) catalog. Once that is suc- 
cessfully done, the second pass computes the WCS rel- 
ative to a source catalog made from a target exposure. 
All exposures of the same target are fit relative to the 
same reference, regardless of the filter used. Generally 
for this study, the stars used in the final WCS solutions 
number in the two- to three-thousands, which provides 
astrometry accurate typically to Of! 02 rms. 

The final product of the SAO Megacam reduction 
pipeline is the multi-extension file for each exposure. We 



36 SExtractor and Swarp are hosted at |http: //www. a stromatic . 
|net/| 
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run Swarp to mosaic and coadd exposures for each tar- 
get and filter, weighting by the relative zeropoints of each 
frame. 

3.3. Photometry 

SExtractor is used to find sources in the images and 
perform photometry. We operate SExtractor in dual- 
image mode, with the r'-band image serving as the de- 
tection image. In this mode, objects are detected in the 
r'-band image while photometry is done in the g' , r', 
or i'-band image, and these final catalogs are joined to 
produce a catalog of colors and magnitudes in the in- 
strumental system. We use MAG_AUT0 photometry. The 
colors of stars and galaxies a re then calibrated using Stel- 
lar Locus Regression (SLR, Hig h et aL||2009[ ). SLR cali- 
brates colors by fitting the instrumental stellar locus to 
that of ~ 10 5 stars in the SDSS. Cross matching with the 
2MASS allows us to solve for the zeropoints of individ- 
ual bands as well to produce calibrated magnitudes. The 
resulting photometry is effectively dereddened, as this is 
an inherent feature of the method, so we do not apply 
any additional Galactic extinction corrections. 

Transforming Clay-Megacam photometry to the SDSS 
system with SLR requires estimating color terms. We 
use gri-b&nd photometry of SPT-CL J0516-5430 on the 
SDSS photometric system, acquired with IMACS on 
Magellan Baade, operated in f/2 ima ging mode. Thes e 
IMACS data are used in the analysis of High ct al.|pOTo]) , 
and are described in detail in that work. We cross-match 
point sources with the IMACS catalogs, obtaining the 
color transformations 



ffCMC - 3SDSS = C g + 0.100(.9 - i)sdss 
fCMC - r S uss =C r - 0.022(# - i) SD ss 
«cmc - *sdss = C t - 0.025(g - Osdss, 



(1) 
(2) 
(3) 

where zeropoints C are nuisance parameters left free in 
the fit in addition to the slopes, and CMC denotes Clay- 
Megacam 37 . These measurements arc illustrated in Fig- 
ure [I] We chose the multiplier g — i because it gives 
maximal leverage on the color-term measurement over 
the full range of stellar temperatures, and thus over the 
catalog's color space. 

The systematic uncertainty of the color calibration is 
estimated at 0.01 mag to 0.03 mag in g — r and r — i, 
likely dominated by the non- uniformity of Galactic ex- 
tinction, flat-field ing, and details of the photometry (see 
High et al. 2009). Magnitude calibrations are uncertain 
at ~ 0.05 mag m all bands, dominated by the overall ac- 
curacy of the 2MASS point-source catalog. These levels 
of uncertainty are sufficient for our purposes. 

Table [2] lists the basic imaging information for the five 
clusters. Depths are estimated from the median magni- 
tude of sources whose signal-to- noise ratio is 5. Seeing is 
estimated in the WL band, r'. 

We test the effect of photometric zeropoint errors on 
the final WL-SZ mass ratios. Systematic errors in pho- 
tometry enter into the mass analysis through the estima- 
tion of the critical surface density for each cluster (Sec- 
tion \4~2l. We estimate this quantity using publicly avail- 



able photometric redshift catalogs of external fields, to 




1.0 1.5 2.0 

(9 " Osdss 
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Fig. 1. — Clay-Megacam color-term measurement. Data points 
are the photometry of point sources that have been matched be- 
tween the Clay-Megacam and IMACS catalogs, the latter of which 
has already been transformed to the SDSS system. The best-fit 
lines, whose slopes are equal to the color-term coefficients, are also 
shown. 



TABLE 2 
Imaging data 



Cluster Name 


Mag 


nitude of 5tr 


Point Source 


r Seeing 




9 


r 


i 




SPT-CL J0516-5430 


27.1 


26.7 


26.0 


0'.'67 


SPT-CL J2022-6323 


26.3 


26.2 


25.3 


0'.'88 


SPT-CL J2030-5638 


26.4 


26.2 


25.4 


0'.'80 


SPT-CL J2032-5627 


25.9 


25.8 


24.4 


0'.'82 


SPT-CL J2135-5726 


26.6 


26.1 


25.4 


0'.'89 



37 We emphasize that gcMC = 
and so forth for the other bands. 



9' and gsDSS = 9 in this work, 



Note. — Bas ic properties of the imaging of the five clusters. 
Column 1: R12 cluster designation. 

Columns 2—4: the median magnitude of sources whose signal-to- 
noise ratio is 5. 

Column 5: the median of the stellar FWHM across the entire 
coadded image. 

which we apply the same photometric cuts as are ap- 
plied to the Megacam catalogs. If there is an offset in 
the Megacam z'-band zeropoint relative to that of the 
standard catalog, then we effectively probe a population 
that is different than that from which we infer the source 
redshift distribution. We test the effect of photomet- 
ric error, 5i, by cutting the standard photo- z catalogs 
at icfht > 24 — Si and repeating the full analysis. The 
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level of photometric accuracy estimated here (5%) causes 
changes in the WL-SZ mass ratios at the sub-percent 
level, which is significantly subdominant to the statisti- 
cal uncertainty and the largest systematic uncertainties. 

4. CREATING SHEAR CATALOGS 

In this section, we summarize the standard theoret- 
ical framework on which WL mass measurements rest, 
including the two primary quantities that must be es- 
timated from data: the critical surface density and the 
reduced shear. 

4.1. Tangential Shear 

Weak gravitational lensing of extended sources by 
spherically symmetric mass overdensities induces a mean 
shear in a direction oriented tangentially to the center of 
mass. Tangential shear, j+, is calculated from the Carte- 
sian components of shear, (71,72), as 

7+ = -7i cos(20) - 72 sin(20) (4) 

(see for example Mellier|1999 Section 2, for an overview 
of WL shear). Indices i <G {1, 2} correspond to horizontal 
and vertical image coordinates. Here, 71 is the compo- 
nent along the horizontal axis (position angle = 0°) 
and 72 is the shear at position angle <fi = 45°. Cross 
shear is calculated as 



7> 



-71 sin(2</>) + 72 cos(2<?!>); 



(5) 



this is the shear component oriented at 45° with respect 
to 7+. The azimuthally averaged cross shear (7 
a function of radius provides a diagnostic for residual 
systematics, because no astrophysical effects, including 
lensing, produce such a signal. As a consequence, a non- 
zero (7x ) indicates the presence of some types of residual 
systematic error, though we note that this is not an ex- 
haustive test. 

The mean tangential shear as a function of radial dis- 
tance in the plane of the sky at the cluster redshift, R, 
depends on the projected surface density, E(i?), as 



<E) (<#)-£(#) 



^crit 



(G) 



( Miralda-Escude|1995 ) . This depends on the critical sur- 
face density, 

c 2 1 

where c is the speed of light, G is the gravitational con- 
stant, and P = D[ s /D s is the lensing efficiency. Quan- 
tities D are angular-diameter distances, and I indicates 
the lens (the cluster) while s indicates sources. 

The observable quantity is not the shear but the re- 
duced shear, g, which relates to the shear as 



7 = (1 - K)g 



(8) 



via the convergence, n = E/E cr i t . Estimating true shear 
therefore requires estimating convergence, whic h we per- 
form in radial bins, as discussed in Section [5.1| 

The two key quantities in measuring mass from WL 
data are thus the critical surface density and the reduced 
shear. We describe these two steps in the following sec- 
tions. 



4.2. Cluster- galaxy Decontamination and Critical 
Surface Density Estimation 

Calculating WL masses requires estimating the criti- 
cal surface density (Equation ([7])), which is a geomet- 
ric quantity containing ratios of angular diameter dis- 
tances between the observer, the lens, and sources. This 
requires redshift information for the lens and sources. 
Three bands are not sufficient for estimating source 
redshifts in the Clay-Megacam data themselves, so we 
rely on photometric redshift catalogs of other, non- 
overlapping surveys that have integrated to equal or 
greater depths, under the assumption that the mean un- 
derlying galaxy population is the same everywhere in the 
sky. The Canada-France-Hawaii Telescope Legacy Sur- 
vey (CFHTLS) Deep field ca talogs are sufficient for this 
purpose ( |Coupon et al. 2009). For the E cr i t estimate from 
the photo- z catalog to accurately reflect the population 
in the Clay-Megacam images on average, cluster galaxies 
must be removed from the Clay-Megacam catalogs. 

We remove cluster galaxies from the shear catalogs us- 
ing the same CFHTLS catalogs as a guide. We first 
plot the density of galaxies at 19 < ?cfht < 25 in the 
CFHTLS-Deep catalogs in (g — r, r — i) color-color space. 
Galaxies with photometric redshifts of |z p hot — zi \ < 0.05, 
i.e., near the cluster redshift, are considered contami- 
nants or "non-sources" , and all other galaxies are consid- 
ered "sources". The densities of these two populations 
are shown in Figure [2] We then define a simple polygon 
that encloses the majority of the non-source galaxies. 

Because the clusters are at different redshifts, the lo- 
cation and shape of this polygon is a function of z\ 
in general. Rather than defining five different poly- 
gons, we only define two: one that excises cluster galax- 
ies for 0.28 < 2; < 0.35 (two clusters), and one for 
0.35 < zi < 0.43 (three clusters). Figure [3] shows the 
photometric redshift distribution after removing galax- 
ies within the high-z; polygon. We test the contamina- 
tion after cuts by using these polygons to remove galax- 
ies from the CFHTLS-Deep catalogs, and measuring the 
fraction of galaxies with photometric redshifts satisfying 
l^phot — Zi I < 0.05. Under this test, the fraction of con- 
taminants is < 0.1% for all cluster redshifts considered 
here. 

After transforming the Clay-Megacam catalogs to the 
CFHT-Megacam photo metric system using color terms 



(Regnault et al. 
and color cuts to 



20091), we apply the same magnitude 
he Clay-Megacam catalogs. The radial 



distribution of catalog galaxies before and after the color 
cuts is shown in Figure [4] for a representative cluster. 
This procedure removes most of the radially decreasing 
trend of galaxy densities such that the final galaxy sur- 
face density is uncorrelated with radius. This constitutes 
some empirical evidence that the procedure is removing 
cluster galaxies. Some residual surface density trend with 
radius may be expected due t o WL magnification (e.g. 
Bartelmann fc Schneider|2001[ ). We measure the slope of 
dlogiV/dm of selected sources to be ~ 0.4 in our data, 
so the effect would be small. 

We estimate (/3) and (/3 2 ) from the photo-z distribu- 
tion in the CFHTLS-Deep catalogs after these selections 
are made (Figure pi) . The f3 distribution is expected to 
vary between fields, both in the CFHTLS Deep fields 
and the SPT cluster fields, due to finite galaxy counts 
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Fig. 2. — Color-color distribution of 19 < icFHT < 25 galaxies in the CFHTLS-Deep photometric redshift catalog. Upper left: the blue 
region denotes the distribution of source galaxies (|.Zphot — £l| > 0.05), and red denotes non-source galaxies (|z p hot — £l| < 0.05), both 
for an example cluster in the high redshift bin, SPT-CL J2022-6323 at z; = 0.383. Contours delineate isodensitics in color— color bins of 
size 0.1 mag X 0.1 mag, and are logarithmically spaced by factors of 10 ' 5 starting at 10. The g ray polygon delineates source/non-source 
regions. The yellow circle shows the typical color of a luminous red galaxy at the cluster redshift | |Lopes|2007| . The remaining panels show 
one-dimensional histograms of projections of these colors. These panels illustrate that one-dimensional color cuts are a less efficient way of 
removing contaminant galaxies. 



and cosmic variance. We estimate the uncertainty this 
induces in the mean mass ratio results by repeating the 
entire analysis using each CFHTLS Deep field individ- 
ually, and averaging the results. The uncertainties on 
resulting WL to SZ mass ratios are < 2%, a level that is 
subdominant to the overall statistical uncertainties and 
other systematics. We adopt the mean values of (/3) and 
(f3 2 ) over the four CFHTLS Deep fields in the baseline 
analysis. These are reported in Table [3j 

We perform consistency checks on the WL mass anal- 
ysis by modifying this photometric selection procedure 
in two ways. First, we test defining color-color poly- 
gons that remove the vast majority of objects at z p hot < 
zi +0.1, i.e., both near and in front of the cluster. This 
significantly widens the area covered by the polygon, and 
the effect is that a large number of galaxies that are 



in fact behind the cluster are eliminated as well. This 
roughly halves n ga \ and does not change (/3) estimates 
significantly. The WL-SZ ratios are consistent with unity 
and with the baseline result under this test. 

Second, we test using limits of zcfht = 24.5 and 24, 
which are brighter than our baseline limit of icfht = 25. 
The effect is to probe a brighter mean population of 
galaxies in the catalogs, which causes (/3) and n ga \ to 
both take smaller values. This increases the statistical 
uncertainty of the WL masses, but restricts the catalogs 
to magnitude regimes in which the CFHTLS-Deep pho- 
tometric redsliif^cj^takogshave been explicitly tested and 
verified ( Coupon et al.|20 09). The resulting WL-SZ aper- 
ture mass ratios in both cases are statistically consistent 
with the baseline result, and with unity. 
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TABLE 3 
Source galaxy properties. 



Cluster Name 




(P) 


</3 2 > 


"gal 










(arcmin -2 ) 


SPT-CL J0516-5430 


0.294(1) 


0.64 


0.43 


14.8 


SPT-CL J2022-6323 


0.383(1) 


0.54 


0.32 


10.6 


SPT-CL J2030-5638 


0.40(4) 


0.53 


0.31 


13.3 


SPT-CL J2032-5627 


0.284(1) 


0.65 


0.44 


9.0 


SPT-CL J2135-5726 


0.427(1) 


0.50 


0.29 


12.2 



Note. — This table summarizes the basic properties of 
the sources used in the weak-lensing mass measure ment after 
making the catalog selections described in Section |4.2| 
Column 1: R12 cluster designation. 

Column 2: cluster redshift. The number in parentheses is 
the uncertainty in the last digit. 



Column 3 
Column 4 
Column 5 



the mean of /3 = Di s /D s after cuts, 
the variance of /3 after cuts, 
the source number density after cuts. 
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Fig. 3.— CFHTLS Deep field 3 photo-z distribution of 19 < 
*CFHT < 25 galaxies after we have identified sources and non- 
sources using the color cuts illustrated in Figure [2] A redshift of 
0.383, corresponding to SPT-CL J2022-6323, is shown. 



4.3. Shear Measurement Method 

The second key ingredient in estimating cluster masses 
in WL analyses is reduced shear. To estimate reduced 
shear in the r'-band images we emplo y the met hod of 



Kaiser, Squires, fc Broadhurst 
and |Luppino &; Kaiser (1 1997b , 



r id [Cu pping & Kaiser ( |1997p , with the modifications 
of |rloekstra et al.| ( |l998 ). This method bases shear mea- 
surements on Gaussian- weighted second-order image mo- 
ments, 



(1995 hereafter KSB) 



d 2 xXiXj W(x)f(x), 



0) 



where W is the Gaussian function, / is the object image, 
and the origin of the coordinate system has been itera- 
tively determined from the first-order moment using the 
same weight. Variables Xi (i £ {1,2}) indicate vertical 
and horizontal image coordinates. The polarization is a 
spin 2 pseudo- vector e a (a £ {1,2}), where 



pi 



hi - h 
hi + h 



and 62 



2/l2 



hi 



(10) 



9 [arcmin] 

Fig. 4. — Radial distribution of galaxies in the Clay-Megacam 
catalog of cluster SPT-CL J202 2-63 23, before and after making the 
color cuts described in Section |4.2| These data are representative 
of all clusters. The color cuts reduce the clustering signal such 
that the final catalog used in the weak lensing analysis is roughly 
flat. WL masses are computed outside the shaded region, from 8i 
(vertical dashed line) to 02 = 12'. 

The magnitude of the polarization is e = (e 2 + e?,) 1 ^ 2 . 

The resulting polarizations must be corrected for the 
effects of the anisotropic smearing by the PSF. To this 
end, we fit fourth order polynomials to the stellar P* 
and Pf™ (the diagonal e ntries of the shear- and smear- 
polarizability tensors, see KSB ) and remove outliers. To 
the surviving objects we tit a fourth order polynomial to 
both ei and e-i. We do this for a range of weight func- 
tions, which is necessary fo r computing the pre- seeing 



shear polarizability P 7 (see Hoekstra et al. 1998). The 
rms of the residuals for each polarization component af- 
ter the fourth order fit are < 0.005 for all clusters. These 
are small compared to the magnitude of signals we are 
seeking, which are about 0.01 to 0.1. While residuals 
of this size can in principle lead to large shear bias for 
any given faint galaxy locally, we expect this effect to 
be zero on average in the radial bins because the spatial 
residual pattern is consistent with noise. A representa- 
tive example PSF polarization plot is shown in Figure 
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[5j This model is used to correct for PSF anisotro py us 



ing the procedure de scribed in the original works (KSB 
|Hoekstra et~aL|l998 l 



1'he next step is to account for the seeing, which lowers 
the observed shear signal. This is done by r escaling the 
polarizations to their "true" value s by P 1 (Luppino & 



Kaiser 1997 Hoekstra et al. 19981. The measurements 
of P' 1 for individual galaxies are too noisy, and instead 
we bin the measurements in galaxy size and use this to 
compute the value as a function of size. The reduced 
shear g is then estimated as 



9i = ei/P^ 



(11) 



This implementation of the KSB| shear analysi s pipeline 
is described in m ore detail in previous works (Hoekstra 
et al.||1998[ |2000|) . It wa s subject to blind tests' by the 



Shear Testi ng Program (Heymans et al. 2006 Massey 



et al. 2007), wherein it achieved accuracy at the O.O27 
level under low PSF anisotropies, such as Clay-Megacam 
exhibits. This level of bias induces mass errors that are 
subdominant to the statistical and systematic uncertain- 
ties. 

5. WEAK-LENSING MASS MEASUREMENT 

In this section we discuss how we estimate mass from 
the shear catalogs. We employ two types of WL mass 
estimators, both of which make use of the same, az- 
imuthally averaged shear profiles. We test the accuracy 
of our algorithms using ray-traced iV-body simulations, 
and we also describe the WL convergence field recon- 
struction procedure. 

5.1. Shear Profiles and Corrections 

We compute binned shear profiles using a weighted av- 
erage in each radial bin. The weight for each galaxy is 



pi- 



1 



(12) 



where a~ t is the scatter in shear due to the intrinsic el- 



Leau- 



lipticity of the ga laxies (set to 0.3; see for example 
thaud et al. | 2007 ) and (Ae 2 ) is the error estimate for the 
polarization (Hoek stra et al.|[2000 |. The weighted mean 

{+-xF 



for i G 



is then 



(7.) 



J2n w nji,n 



and errors on the mean, er( 7 . 

1 



, are computed via 



'(7<> 



(13) 



(14) 



Two corrections are applied to the binned shear data. 
The first correction accounts for a known error in the 



binned shear data due to the averaging operation. Seitz| 
& Schneider ( 1997[ ) show that estimating the critical sur- 
face density of each cluster using the mean of the j3 dis- 
tribution in redshift catalogs induces an error in the ob- 
served reduced-shear, which can be corrected to first or- 
der using 



(<?, 



'obs 



(Str 



1 + 



f(f3 2 



1 K. 



(15) 



We adopt the model for k described below to make this 
correction. 

The second correction is to transform the reduced shear 
g to shear 7 by estimating the n in each bin. We accom- 
plish this by jointly fitting for the shear and the conver- 
gence profiles to those predic ted by Na varro, Frenk, &] 
Whitel (|1997l hereafter |NFW l) den sity profiles assuming 

D08l) 



Dufiy et a. 
redshift sea 



hereafter | 
rerein the I 



concentration— mass- 
erived mass is used as 



the input. The NFW family of density profiles take the 
form 

P{r) = (rc/r 20 o)(l+V C /r 20 o) 2 ' (16) 

where r is the three-dimensional radial distance, r2oo is 
the radius at which the mean NFW overdensity is 200 
times /9 cr it at the cluster redshift, c is concentration, and 



200 



3 ln(l 



c/(l 



(17) 



This fitting procedure result s in e stimates of the 7 pro- 
file under the assumption of pD08 concentrations, and it 
also yields the convergence profile used to perform the 
redshift distribution correction in Equation [T5] as well as 
total spherical masses in our main results. 

Observed shear profiles, after all of these corrections 
are applied, are presented in the Appendix. The cross 
shear is consistent with zero, as expected for a signal 
that is not appreciably contaminated by systematics. 

5.2. Aperture Masses 

One of the advantages of the tangential shear is that 
it can be used to directly constrain the projected mass 
within an aperture. Aperture masses are integrals 
over the shear data multip lied by a filter (see Fahlman 
et al.||1994| |Schneider||1996| . These aperture masses are 
parameter-free in the sense that the observable quantity 
gives cylindrical mass constraints without reference to 
analytic density profiles. If the filter is compensated, 
they are also insensitive to the mass-sheet degeneracy. 
The classical ^-statistic for circular apertures, which uses 
a particular choice of filter, is 38 



C(0i) 



KOi 



Wfli<e<e 2 



1 oi/oi 



de 



<7+>(f) 



(18) 



( Fahlman et al.|1994| Kaiser 1995). The aperture mass 
is then 

M ap (0x) = 7r(^A) 2 S crit C(^i). (19) 

The square of the measurement uncertainty of the aper- 
ture mass statistic takes an analytic form, 



1 -oi/oi 



AO 



(7+) 
6» 2 



(20) 



is given in Equation (14). In addition 



where er( 7+ ) 

to the measurement uncertainty, LSS along the line of 
sight contributes a random uncertainty to the WL aper- 
ture masses, which we add to the formal statistical WL 

38 This £ statistic is in general a function of two angles, 9i and 
62, but we fix 82 for all analyses and so have simplified the notation. 
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Measured 



Residual w.r.t. fit 



12 
10 



Q. 
O 



5% polarization: 




2 4 6 8 10 12 
a; [10 3 pixels] 

Fig. 5. — PSF polarization plot of the coadded r-band image of SPT-CL J2022-6323, which is representative of the sample. Left: 
polarization magnitude and direction, as measured from the Gaussian-weighted image moments of point sources. Right: same as the left 
panel, after a fourth order polynomial in the two spatial dimensions is fit to each component of polarization. Residuals between the data 
and the model are < 0.005 rms for all clusters. 



m ass uncertainties in qu adrature using the prescription 



of Hoekstra et al. ( 2011[ ). LSS uncertainties are 15% to 
20% tor these clusters. 

We fix 6 2 to 12' for all analyses, which is roughly the 
maximum radius of the Megacam imaging. We set 0\ to 
^?500,sz/ Di- 

Aperture masses provide lower limits on total cylindri- 
cal masses within 0±. In the limit of 02 — > oo, the aper- 
ture mass converges to true integrated cylindrical mass 
within 8i , but real data extend to finite radius, so 02 must 
be finite. Direct comparisons of the observed aperture 
mass to other spherical or cylindrical mass-observables 
therefore requires care. 

5.3. Spherical- to Aperture-mass Transformations 

Some extra computation makes direct comparisons be- 
tween aperture masses and spherical masses possible. 
This is accomplished by assuming an analytic three- 
dimensional profile consistent with a given spherical mass 
estimate, projecting it to two dimensions, calculating the 
shear assuming some E C rit, and then computing the aper- 
ture mass statistic given this predicted shear. 

To illustrate this procedure concretely, say we are 
given some spherical mass estimate, MA, bs, derived 
from another observable such as the SZ effect. We 
wish to compare this to the WL aperture mass at 
the same radius, M aPi wL(-RA,obs), where i?A,obs = 
[3M"A,obs/ (A 47T/9 cr it (^))1 1 ^ 3 . We first assume a three- 
dimensional |NFW| profile, with total spherical mass 
within i?A,obs equal to MA >0 bs, and concentration c taken 
from the mean of clusters of these masses and redshifts 



dius is 0i = i?A,obs/A- The result, M aP:0hs (R A , ohs ), 
is the aperture-equivalent mass of AfA.obs, which was 
determined by some other method or observable. It is 
the direct analog of the WL aperture mass measured at 
the same radius, M aPi wL(-RA,obs), such that the ratio of 
the WL to this aperture-equivalent quantity is unity in 
the absence of other systematic errors. We perform this 
transformation on SZ mass estimates to test for statisti- 
cal consistency with the WL aperture masses. 

We propagate spherical mass uncertainties to aper- 
ture mass equivalent by computing the aperture mass 
of MA,obs ± &m numerically using the same proce- 
dure. As illustration, aperture mass uncertainties are 
(0.54, 0.67, 0.81) times the values of their spherical 



NFW counterparts for c = (1,3,10) at these clusters' 
typical -R500 radii. 

5.4. Calibration Tests with Mock Catalogs 

A number of recent works have shown using JV-body 
simulations that WL derived mass estimates are biased 
at roughly the —5% to —10% level (|Becker fc Kravtsov 



2011 Bahe et al. 2012 



Rasia et aT 
of the 



WT 



"We" 



per- 



according to the D08 c-M-z scaling relation. The corre- 
sponding two-dimensional projection o f this profile and 
predicted shear take an analytic form (Wright & Brain- 
erd 2000). We use the same £ cr jt as used in the WL 
analysis. Then, from this shear we compute (, again 
using the same filter as used on the WL data. Specif- 
ically, the outer radius, 02 = 12', has been set by the 
size of the Megacam imaging footprint, and the inner ra- 



torm similar calibration tests of the WL mass statis- 
tics presented in our work using mock shear catalogs of 
220 deg 2 of sky. The catalogs are drawn from an iV-body 
dark matter simulation of a standard ACDM universe. 
The dark matter halos are populated with galaxies using 
ADDGALS such that they repro duce known lum inos- 
ity, color, and clustering relations (Wechsler 2004, R. H. 
Wechsler et al., in preparation). Shears are assigned to 
each galaxy by raytracing through the A-body simula- 
tion 39 . 

We compute the masses of 280 A-body halos at red- 
shifts 0.25 < z < 0.45 with masses M 200 > 10 14 h^M Q . 
To reflect the choices we have made in analyzing the 
Clay-Megacam data, we only use shear profiles between 

39 The simulated shear catalogs were kindly made available to 
use by R. Wechsler, M. Busha, and M. R. Becker. 
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i?50o and 12', where i?soo is extrapolated using the |D08 
concentration-mass-redshift scaling relation with M200 
determined from the TV-body halo finder as the input. In 
order to minimize statistical uncertainty in these tests, 
we use perfect knowledge of the shear field at the galaxy 
locations (i.e., we do not include intrinsic shape noise) 
as well as perfect source redshifts, and we select source 
galaxies as those at \z — z\\ > 0.05. The resulting bias 
between measured masses and appropriately transformed 
A^-body masses, averaged over the sample, ranges from 
—6% to —13%, consistent in magnitude and sign with 
the previous works. These tests carry a statistical uncer- 
tainty of abou t 2%. The precise bias values are reported 
in Section pO} We do not apply these bias corrections to 



any mass estimates presented in this work. 

5.5. Convergence Field Reconstruction 
We reconstruct the convergence field i n two dimensions 



using the method of|Kaiser & Squires| (|1993J). The pur- 

ik of the re- 



pose is to test the impact ot using the peak 
constructed k field as the cluster center when computing 
shear profiles. The convergence is estimated up to a con- 
stant as a sum over galaxies n, 



k(x ) = -- 



1 



"galTT 



y^XaO?n - x )2e a /P 1 



(21) 



where n ga i is the mean surface density of source galaxies 
and Xa = {x 2 - y 2 ,2xy}/9 4 and 6 = {x 2 + y 2 ) 1 / 2 . The 
constant is set such that the mean convergence across 
the full field is zero. The wide field-of-view of the ob- 
servations allows us to avoid artifacts from this method 
which have caused problems in the past. 

We pixelate the convergence maps at O^pixeP 1 . The 
errors per pixel are independent, and are computed 
as (Jpix ~ <t k / -y/n ga i0.'4 2 , where a K ~ er 7 ~ 0.3 and 
n ga i is measured from the data and assigned units of 
arcmin" 2 . We then smooth the map with a Gaussian 
of size Csmooth = 3'. The peak in the smoothed con- 
vergence field, max(K), is identified as the maximum 
value across the entire smoothed map. The signal- 
to-noise ratio in the smoothed peak value is S/N = 
max(K)/(<7pi X / •y/27TCT smoo th/0.'4). The uncertainty on 
the peak position is then FWHM slnoot i,/(S/JV), where 

FWHM smooth = v / 81n2a smooth - This is equal to - 45" 
for all five clusters. Using the convergence field peaks as 
the cluster centers gives a mean ratio of WL-SZ masses 
that is consistent with unity and with the baseline result 
to within the statistical uncertainty. 

6. RESULTS 

We report WL-derived masses, then te st th e overall 
accuracy of the SZ mass determination of |R12| by mea- 
suring the mean ratio of equivalent WL and SZ mass 
estimators. In the baseline analysis we 

• select source galaxies at 19 < icfht < 25 with 
colors exterior to the polygon shown in Section 4.2 



use concentrations of |D08[ and 

use SZ centroid as the cluster centers. 



We then alter various steps in the analysis in order to 
test the robustness of the result, as well as to estimate 
the magnitude of various potential sources of error. Shear 
profiles and aperture mass profiles are presented in the 
Appendix, in addition to optical, SZ, and convergence 
maps. 

6.1. Masses and Mass Ratios 

In this section we use three different methods for esti- 
mating mass from these data. In all cases, we only use 
shear data at i?5oo,Sz/-D 1 < < 12'. 

In the first method, we measure aperture masses with 
9i = i?5oo,Sz/A and 62 = 12'. Following the procedure 



outlined in Section 5.3 we derive the equivalent estima- 



tor from the SZ masses by assigning an NFW profile con- 
sistent with M5oo,sz, computing the predicted shear for 
such a profile, and then calculating the aperture mass 
from the predicted profile over the same radii. These 
masses are inherently projected, two-dimensional quan- 
tities. We plot the results in the top panel of Figure 
1 

In the second method, we estimate the spherical 
mass by fitting NFW profiles to binned shear data 
at ^500,sz/A < < 12' and computing the total 
mass of the best fit profile at i?5oo,sz- We label this 
-Mnfw,wl(-R500,Sz)- The results are to be compared to 
M5oo,sz directly. Because of the use of a radius from an 
external source, the radius at which the mass is quoted 
is not the radius where A = 500 (the overdensity factor 
with respect to the cosmological critical density) in the 
best-fit NFW model. These masses are plotted in the 
bottom panel of Figure [6] 

And in the third method, we use the same best-fit 
profiles above to estimate the spherical mass within the 
radius where A = 500 as determined from the best-fit 
model itself, i.e., M 5 oo,wl- This is also to be compared 
to M500 sz directly. These masses are plotted in Figure 

m 

The expectation of the ratio of WL to SZ masses is 
unity for all of these methods in the absence of sys- 
tematic errors from, for example, the concentration and 
cluster centering assumptions. We note that we use a 
concentration-mass-redshift scaling relation in all ratio 
tests: for the aperture mass comparison, it is only used 
to transform the SZ-derived mass to an aperture mass 
equivalent, whereas for both spherical mass comparisons 
it is only used in the WL shear profile fit. 

The mass results and derived quantities are reported 
in Table |U 

The mean ratios of the three WL to SZ mass statis- 
tics are summarized in Table [5] In all cases we report 
weighted mean values, where, for each cluster n, weights 
l/cr„ are a combination of the WL aperture mass statis- 
tical uncertainty (including the estimated LSS contribu- 
tion, which is between 15% and 20% fo r the se clusters) 
and the total SZ mass uncertainties from |R1 2] propagated 
to the derived quantity when necessary. The uncertainty 
on the mean is computed via 1/er 2 = J~) n 1/c 2 . 

We note that this method does not take into account 
any correlated uncertainty between clusters. The SZ 
mass estimates from R12 have been marginalized over 
cosmological and scaling relation parameters, and this 
results in a ~ 10% systematic uncertainty that is highly 
correlated between clusters. We have checked for the ei- 
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Fig. 6. — Top: weak-lensing aperture masses vs. inferred aperture masses from Sunyaev-Zel'dovich effec t data. The one-to-one line is 
also shown. Bottom: spherical WL masses vs spherical SZ masses. The WL masses are determined by fittingjNFW profiles to weak lensing 
shear data at -R500 Sz/Di < 9 < 12' and evaluating the resulting NFW mass profile at -R500 SZ- 



feet of these correlations. Briefly, we adopt a Gaussian 
likelihood for both the SZ and WL mass statistics and in- 
troduce nuisance parameters representing the true mass 
of each cluster. We introduce an additional free param- 
eter representing an overall scaling of the SZ mass esti- 
mates. We then explore the resulting likelihood surface 
using an MCMC. If we use a diagonal covariance matrix 
with the uncertainties given in Table |4j we recover the 
nominal mean ratios of WL to SZ mass statistics of Ta- 
ble [5] at the maximum likelihood points in the chains. If 



we use the covariance matrix for the five clusters as esti- 



mated from the R12 cosmological chains, we find that the 
increase in the uncertainty on the mean ratios is small 
(~ 10% of the baseline ratio uncertainty values) when 
compared to when we s et these correlations to zero. As 
discussed in Sect ion [(O] WL systematic uncertainties are 
estimated to be significantly smaller than the statistical 
components, so we only use a diagonal WL covariance 
matrix. Thus, we ignore the correlated component to 
the SZ and WL cluster mass estimate uncertainties for 
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TABLE 4 

Cluster mass results 



Cluster Name 


^500, SZ 


A^ap,WL(-f?500,Sz) 


A^ap,Sz(-f?500,Sz) 


A^NFW.WL (#500, SZ ) 


A^500,WL 


*^500,SZ 




(/i _1 Mpc) 


(lO 14 /!"^©) 


(lo 14 ^- 1 ^®) 


(lO^/l-iMQ) 


(lo^/i-iJUe) 


(10 14 /i- 1 Af Q ) 


SPT-CL J0516-5430 


0.84 


4.70 ± 1.49 


5.31 ± 0.60 


4.43 ± 1.16 


4.35 ± 1.57 


4.60 ±0.90 


SPT-CL J2022-6323 


0.69 


3.03 ± 1.64 


3.42 ± 0.41 


3.08 ± 1.12 


3.27 ± 1.64 


2.75 ±0.60 


SPT-CL J2030-5638 


0.64 


2.62 ± 1.43 


2.87 ± 0.39 


1.93 ±0.82 


1.77 ± 1.15 


2.28 ±0.57 


SPT-CL J2032-5627 


0.77 


6.62 ± 1.66 


3.99 ± 0.47 


5.75 ± 1.43 


7.71 ± 2.23 


3.40 ± 0.70 


SPT-CL J2135-5726 


0.77 


5.21 ± 1.88 


4.99 ± 0.52 


4.37 ± 1.29 


4.57 ± 1.87 


4.03 ± 0.78 



Note. — Results from the WL measurements, as well SZ data and derived quantities. We have divided two dimensional aperture 
masses fro m sph erical masses. 
Column 1: |R12| cluster designation. 

Column 2: cluster radius determined from the SPT SZ mass estimate in Column 7. 
Column 3: WL aperture mass at -Rsoo.SZ- 
Column 4: SZ aperture mass at Rsoo.SZ- 

Column 5: WL mass at i?5oo,SZ from | NFW| p rofile fitted to shear data at i?soo,Sz/f ( < 9 < 12'. 
Column 6: WL mass at -Rsoo.WL f rom IN Fw I profile fitted to shear data at i?5oo,Sz/£* I < 9 < 12'. 
Column 7: SZ mass estimate from IRlzl 



TABLE 5 

Mass ratio results and consistency tests 



Procedure 


Aperture 


Masses 


Spherical Masses at -Rsoo.SZ 


Spherical Masses at -Rsoo.WL 




Mean Ratio 


Total Scatter 


Mean Ratio 


Total Scatter 


Mean Ratio 


Total Scatter 






(96) 




[%] 




(%) 


Baseline results 


1.04 ±0.18 


33 


1.07 ±0.18 


33 


1.10 ±0.24 


58 


«CFHT < 24.5 


1.05 ±0.20 


37 


1.05 ±0.19 


37 


1.08 ±0.25 


65 


«CFHT < 24.0 


1.07 ±0.23 


31 


1.07 ±0.21 


35 


1.12 ±0.29 


60 


Conservative color cuts 


1.07 ± 0.19 


35 


1.12 ±0.18 


35 


1.16 ±0.23 


58 



Note. — WL-SZ mas s ratio results. The baseline analysis employs a magnitude limit of icFHT = 25 and color cuts 
described in Section 14.21 

Column 1: procedure used to compute mean ratio. 

Column 2: weighted mean of the ratio of WL to SZ aperture masses, (M aPi wL(-R500,Sz)/A^ap,Sz(K5oo,Sz)}. 
Column 3: the unweighted standard deviation of the WL to SZ aperture mass ratio. 

Column 4: weighted mean of the ratio of WL spherical mass (evaluated at -R500 sz) to SZ spherical mass, 
(MjvfW,Wl(-R500,Sz)/A^500,Sz)- 

Column 5: the unweighted standard deviation of the Column 4 ratio statistic. 

Column 6: weighted mean of the ratio of WL spherical mass (evaluated at -Rsoo.wl) to SZ spherical mass, 
(A^500,Wl/M500,SZ>- 

Column 7: the unweighted standard deviation of the Column 6 ratio statistic. 



the remainder of this work. 

The total scatter reported in Table [5] is the unweighted 
standard deviation of the respective ratio data. We do 
not estimate intrinsic scatter because this requires an 
estimate of the level of correlation in intrinsic scatter 
between the WL and SZ mass-observables. Estimating 
these q uant ities is beyond the scope of this work (see 
Section 7.1 for further discussion). 

We have performed additional consistency tests by us- 
ing brighter ?cfht magnitude limits, and by adopting 
more conservative color cuts. The brighter magnitude 
limits are meant to probe the CFHTLS-Deep catalogs at 
magnitudes where the p hoto- z accuracy has been explic- 
itly verified and stated ( Coupon et al. 2009 ). Results of 
these tests are also reported m Table [5} For all tests we 
have performed, mean ratios are statistically consistent 
with unity and with the baseline result to within lcr. 

6.2. Systematic Error Analysis 

We have explored various potential sources of system- 
atic error, as discussed throughout the text, and we re- 
port the effects on the mean WL-SZ mass ratios in Table 
[6] for the dominant or most interesting cases. Changes in 
mean ratios are quoted relative to the baseline results of 
the aperture mass and spherical mass comparisons (Ta- 



ble]!]). 

Our first tests are to measure the effect that the as- 
sumed concentration scaling relation has on the mean 
ratios. Our baseline analysis uses that of |D08| which 
gi ves 3.0 < c < 3.2 fo r these five clusters. The work 
of iMaccio et al.| d2008| ) give 4. 2 < c < 4.5 and lNeto 



etaT. 



Prada et al. (2011 iTmd" 



fl2007| ) give 3~3"<~c < 4.1 

concentration values higher than these other works, with 
50th percentile values of c « 5.0 and 90th percentile val- 
ues of c 7.1. The latter represents the most extreme 
conc entrations we might expect for this sample, while the 
|D08| gives the lowest expected values. This range of con- 
centrations change the aperture mass ratios by 2% to 6%. 
Spherical mass ratios using NFW fits to WL data evalu- 
ated at i?5oo,sz change by 3% to 9%, and those evaluated 
at J?5oo,wl change by 4% to 11%. This level of uncer- 
tainty is less than our statistical uncertainty on the mass 
ratio for the five clusters, but it is the dominant source 
of systematic uncertainty in this work. 

While we use the SZ positions as the cluster centers 
in the baseline analyses, we also test using two other 
definitions of the cluster center: (1) the BCG and (2) 
the peak in the WL-reconstruc t ed n field. We use the 



BCGs selected of Song et al. (2012), which are typi 
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TABLE 6 

Possible sources of systematic error. 



Perturbation 



Change in Mean Ratio 



Aperture Masses Spherical Masses at -Rsoo.SZ Spherical Masses at R$ 



Concentration assumption 



=5 



Duffy et al.N2008> ( baseline) 



iVlacciO ct al. (20 
IMeto et aT7| |Ju(JV|) 
irada et al.fi ^Ull i 50th percentile 
1 2TJTT i 90th percentile 



f rada ct al. 



+0.00 
+0.03 
+0.02 
+0.04 
+0.06 



+0.00 
+0.04 
+0.03 
+0.06 
+0.09 



+0.00 
+0.06 
+0.04 
+0.08 
+0.11 



Cluster center assumption 



Use BCG centers 
Use k map centers 



+0.00 
-0.14 



-0.01 
-0.12 



-0.02 
-0.17 



Other 



Calibration to TV-body simulation 



-0.06 



+0.09 



+0.13 



Note. — Exploration of possible sources of uncertainty in the WL-SZ mass ratios with respect to the baseline results 
shown in Table[5] 



Column 1 
Column 2 
Column 3 
Column 4 



brief description of the perturbation performed, 
change in weighted mean, (M aPiWL (R5oo,sz)/M aPi sz(-R500,Sz)>- 
change in weighted mean, (Af NFWiWL (iJ 5 oo,sz)/A4"500,SZ>- 
change in weighted mean, (Af 5 oo,WL/W 5 oo,sz> ■ 




[ft" 1 M s ] 



Fig. 7. — Spherical WL masses v s spher ical SZ masses. The WL 
masses are determined by fitting NFW profiles to weak lensing 
shear data at R^oosz/Di < 6 < 12' and evaluating the result- 
ing NFW mass profile at iJsoo,WL as determined from the best-fit 
profile itself. 

cally < 30" from the SZ centroid, with the exception 
of SPT-CL J2032-5627, which is nearly 2' away. Mean 
ratio results using the BCG as the cluster center are in 
near perfect agreement with those using the SZ positions. 
However, adoption of the k peak location as the cluster 
center reduces the mean ratios of aperture masses and 
spherical masses by significantly greater amounts. The k 
peak locations have statistical uncertainties of ~ 30" to 
60", which are the largest uncertainties of the three cen- 
ter choices we have considered. While the mean ratios 



using k peak locations are statistically consistent with 
unity and with the baseline results, this represents the 
largest deviation from the baseline. This is probably be- 
cause of noise in centering due to the large statistical 
uncertainties in determining convergence field peak loca- 
tions; indeed, the universal reduction in the ratio values 
using these centers provides some evidence that the con- 
vergence peaks are poorer estimators of cluster centers 
than the SZ and BCG positions given the data, as cen- 
tering errors suppress shear profi les, particularly in the 
inner regions ( |George et al.|[2012| . 



Tests using iV-body simulations (Section [53 result in 
the calibration bias levels reported in Table [5] The sta- 
tistical uncertainties of these calibration bias estimates 
are each 2%. We have not applied these bias corrections 
to any individual mass estimates nor to any mean mass 
ratio statistics presented in this work. We note that such 
bias corrections increase tension with the mean mass ra- 
tio expectation of unity by up to la at most. 

The effects of other known sources of systematic un- 
certainty are smaller than those of concentration. The 
shear measurement method is estimated to be accurate 
to 0.027 in shear (Section]!]), which translates to < 2% in 
the mean aperture mass ratio and < 3% in mean spher- 
ical mass ratios. 
The source redshift distribution inferred from the four 



different CFHTLS Deep fields (Section 4.2) is expected 



to vary due to finite galaxy counts and cosmic variance. 
The source redshift distribution in the SPT cluster fields 
we have observed is also expected to vary across fields in 
the same way. We estimate the uncertainty this induces 
in the mean ratio results by repeating the entire analysis 
using each CFHTLS Deep field individually and aver- 
aging the results. Under the assumption of randomness, 
the mean aperture mass ratio carries an uncertainty from 
finite galaxy counts and cosmic variance of ±1.1%, while 
the mean spherical mass ratios are affected at ±1.0% 
(using i?5oo,Sz) and ±1.6% (using -R 5 oo,wl)- This un- 
certainty can in principle be reduced further by using 
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standard photo- z catalogs of additional, disparate fields, 
or by measuring photo-z's directly in the SPT cluster 
fields, but this is not necessary given the data, as it is 
significantly subdominant to other uncertainties. 

Bias and uncertainty from assumed cosmological pa- 
rameters are also subdominant. We have assumed cos- 
mological paramet er va lues that are slightly different 
from the results of R12 (combining the 100 high purity 
SPT detected clusters with the CMB+BAO+ffo+SNe 
data). To estimate the bias and scatter this in- 
duces, we have test ed m arginalizing the WL to SZ 
mass ratios over the IR12I ACDM chains that combine 
CMB+BAO+iJo+SNe+SPTcL data. The results are 
that both the bias and statistical uncertainty due to 
cosmology are sub-percent in the mass ratio statistics, 
and are therefore significantly below those expected from 
other sources of potential error. 

In summary, no potential sources of systematic error 
that we have explored contributes at or above the level of 
statistical accuracy for this sample of clusters, suggesting 
they are all subdominant. Uncertainty in the assumed 
NFW profile concentration and calibration biases from 
tests using A-body simulations are likely the dominant 
systematic uncertainties. 

7. DISCUSSION 

In this section we discuss the level and possible origin 
of scatter in WL and SZ mass statistics, as well as in- 
dividual cluster systems that exhibit interesting features 
in the data. 

7.1. The Effects of Choice of Radius on Mass Ratios 

We have tested three methods of comparing WL- and 
SZ-derived masses, and find that WL mass statistics eval- 
uated at a fixed radius of i?5oo,SZ exhibit less total scatter 
than spherical WL masses evaluated at i?5oo,WL (Table 

t. This is primarily because i?5oo,WL, effectively de- 
rmined from the best-fit profiles to the WL data, has 
a greater uncertainty than i?5oo.sz- Fixing the radius to 
some externally specified value that has greater precision 
has the benefit of reducing this source of uncertainty. 

On the other hand, masses measured at externally 
specified radii that are themselves determined from a 
mass proxy (such as -ffsoo.sz) must suffer from correlated 
scatter in addition to any that is intrinsic between th e ob- 
servables. For illustration, if the Msoo.sz estimate of R12 
is 20% larger than the true M 500 , then the -Rsoo.sz esti- 
mate scatters up by 7%, and as a result, the WL aperture 
mass and NFW mass evaluated at i?5oo,sz also scatter up 
by some amount, under the assumption of monotonically 
increasing mass profiles. This means that scatter in SZ- 
derived masses induces additional, positively correlated 
scatter between WL and SZ masses due to the use of 
i? 5 oo.sz: in addition to any intrinsic correlation. 

This motivates fixing the radius to some value in arc- 
minutes or in Mpc, or more generally to a value de- 
rived from any data that are minimally correlated with 
the mass-observables or have negligible uncertainty. We 
have tested measuring aperture masses setting 9\ = 
0.8h~ 1 Mpc/Di and 9 2 = 12' for all clusters, where the 
inner radius is roughly equal to the median i?5oo,sz 
for the ensemble. The equivalent statistic for the SZ- 
derived mass is also computed. The weighted mean ratio 
result is (M a p, WL (0.8/i- 1 Mpc)/M ap ,s Z (0.8/i- 1 Mpc)) = 



1.02 ± 0.18, with total scatter of 35%. Compared to the 
baseline spherical mass ratio with 34% scatter in Table [5] 
this increased scatter could be consistent with a ~ 10% 
correlated scatter induced by estimating both the WL 
and SZ spherical masses at i?5oo,sz- 

We have also tested fitting NFW profiles to shear 
data at 0.8/i _1 Mpc/A < 9 < 12', and eval- 
uating the total mass at 0.8ft. _1 Mpc. To com- 
pute the equivalent SZ mass statistic, we evalu- 
ate the NFW profile that is consistent with M 50 o,sz 
also at 0.8ft. _1 Mpc. The weighted mean result 
is (Af NF w,wL(0.8/i- 1 Mpc)/M NFW ,sz(0.8/i- 1 Mpc)) = 
1.12 ± 0.19, with total scatter of 36%. The spherical 
mass ratio result at i?5oo,sz in Table [5] exhibits 15% less 
total scatter (in quadrature) , which again may be due in 
part to the added correlation induced by using i?5oo,sz- 

This preliminary analysis is suggestive that minimum 
variance mass comparisons may be achievable by fixing 
the radius of analysis in Mpc, or some other quantity 
that has negligible uncertainty or is uncorrelated with 
mass. We note that these choice of radii would have 
non-trivial complications for cluster cosmological analy- 
ses, which usually compare the measured cluster abun- 
dance to predictions where cluster masses are typically 
defined as Ma. However, a rigorous comparison of differ- 
ent mass measures and identification of minimum vari- 
ance estimators requires quantifying correlations in scat- 
ter. We leave this to future work. 

7.2. Disturbed Systems and Line-of-sight Structure 

The primary goal of this work is to provide cosmolog- 
ically unbiased tests of the scaling of the SZ o bservable 
with total mass. As described in Section |5.4| we have 
performed calibration tests on mock catalogs based on 
ray-traced A-body simulations for all halos above a given 
mass, as identified by the A-body cluster finder, with- 
out regard to whether the halos have undergone recent 
merger activity or contain significant structures along the 
line of sight. This approximately mimics the SZ selection, 
which is roughly uniform in mass at z > 0.3. The verifica- 
tion exercise indicates no signs of bias to ~ 2% for both 
the aperture mass and spherical mass ratio tests, un- 
der the simple assumptions adopted in that verification 
study. Similarly, we must include any SPT-dctected clus- 
ters that may exhibit merger activity or contain known 
structures in the line of sight, as cutting them out risks 
introducing additional bias in the mean mass ratio tests. 

Nonetheless, identifying disturbed systems and signif- 
icant line of sight structures is of some interest. For ex- 
ample, we note that the inner regions of the shear profiles 
(< Rsoo.sz) of SPT-CL J0516-5430 and SPT-CL J2032- 
5627 show disagreement with the NFW profiles fitted at 
Rboo,Si/Di < 9 < 12', and these two cluster systems also 
show the greatest disagreement between SZ significance 
centroid, BCG, and k peak locations. We explore pos- 
sible explanations here; however, ultimately we conduct 
WL tests in this work in the same way we have tested 
our procedure in A-body simulations, and we therefore 
do not alter the WL-SZ mass comparison methodology 
based on the results of these explorations. 

7.2.1. SPT-CL J0516-5430 

Andersson et al.|(|2011|) note that SPT-CL J05 16-5430 



exhibits north-south elongation in X-ray data obtained 
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with Chandra. Zhang et al. ((2008) also note ellipticity 
in XMM-Newton X-ray data lor this cluster. Both the 
distribution of red cluster-galaxies and the convergence- 
map morphology show north-south elongation as well, 
in rough agreement with the X-ray emission structure 
and SPT SZ significance map. In addition, the conver- 
gence peak disagrees with the SZ significance centroid 
and the BCG locations at about 3er. This suggests that 
this cluster is unrelaxed. It has been shown that core- 
excised X-ray mass-observables are nearly universal for 
relaxed and unrelaxed clusters (e.g., |Mantz et al.||2010a 



Vikhlinin et~al~1|2009a[ ). The SPT SZ observable is simi- 



larly insensitive to the details of the core gas activity, so 
the SZ-derived mass should not be greatly affected. We 
note that the WL- and SZ-derived masses are in excellent 
agreement, which may be due to our use of shear profile 
d ata only outside the c luster core. 

|McInnes et ah] ( |2009[ ) measured the WL mass of SPT- 
CL J0516-5430 using data from the Blanco Cosmology 
Survey. The measured masses agree at the ~ la level 
after analytically accounting for the different concentra- 
tion values used, as well as for the different overdcnsity 
definitions. 

7.2.2. SPT-CL J2032-5627 

The data of the SPT-CL J2032-5627 field show a num- 
ber of interesting features. First, the two brightest clus- 
ter galaxies (both with spectroscopic redshifts consistent 
with 0.28) have relatively large offsets from the SZ cen- 
troid at ~ 1' each, are of similar luminosity, and are 
also consistent with two peaks seen in the WL conver- 
gence map (see Appendix |A|) , which suggests possible 
merging activity. Second, the shear profile is suppressed 
within i? 500 ,sz, similar to SPT-CL J0516-5430. And 
third, there is evidence in the literature for a structure 
in the foreground of the higher redshift, spectroscopi- 
cally confirmed SZ cluster. In this section we explore 
the possible implications of a foreground interloper and 
a plane-of-sky merger on the data. 
The evidence for a foreground cluster or g r oup i s 



strong. As discussed in detail in |Song et al. (2012) 



SPT-CL J2032-5627 is listed by \H.V1\ as having an SZ 
location that is consistent with the optically identified 
cluster A3685 (location 115" ± 150" from the SPT SZ 
position) and the X-ray detected cluster RXC J2032.1- 
5627 (location 87" from the SPT SZ position) A3685 
is assigned z = 0.062 from just one galaxy ( |Fetisova| 
1981| JStruble fc Rood|[T999| , but it is more likely to be 
z = 0.14 because five galaxies near its location are iden- 
ti fied at this redshift in the spectroscopic galaxy catalog 
While RXC J2032. 1-5627 is iden- 



of Guzzo et al 
tified at z 



(|2009 1 
4 from thelGuzzo et al 



(£009} catalog, 
that catalog also contains six galaxies consistent with the 
redshift we identify for SPT-CL J2032-5627, which is de- 
rived f rom 3 1 galaxies from our own observations. |Song| 



et al. 



signa 



(2012) conclude that both the X-ray and the SZ 
s are likely predominantly arising from a massive 
z = 0.28 cluster, but there is also a cluster at z = 0.14 
(possibly the object identified originally as A3685) spa- 
tially consistent with the SZ detection. 

We explore the possible effects of a foreground in- 
terloper and a plane-of-sky merger on the SZ data 
by simulatin g halos using analytic profiles of |Arnaud| 
et al.| ( |2010| ), and on the WL data by simulating NFW 



halos. In our first set of simulations, we model an 
M 2Q0 = W^h^Mo (M 500 = 0.75 x lO^/i^M©) cluster 
at z = 0.14 and an M 500 = 4.30 x lO^h^M® clus- 
ter at z — 0.28. The chosen mass of the interloper 
is motivat ed by the de s ignati on of Richness Class to 
A3685 by jAbell et~al] ( |1989| ), which indicates identifi- 
cation of 30-49 clustered galaxies and is the smallest 
possible richness on that scale. The velocity disper- 
sion of the five galaxies at this redshift in the |Guzzo| 



et al. (2009) catalog (470± 190km s x ) suggests a mass of 

1 M@ using the scaling relation 
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(0. 
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of Evrard et al. ( 2008 ) , which is also consistent. To reflect 
the spatial offset of the convergence peaks seen in the 
data, we separate the cluster centers by ~ 4' by center- 
ing the clusters at the two n peak locations. We confirm 
from these simulations that a foreground object would 
boost WL masses by roughly 50%, while the SZ observ- 
able would be boosted by only < 2%. This configuration 
also reproduces the double peak in the convergence field 
and the suppression of the shear profiles within i?5oo- 
The SZ centroid in this case would be centered on the 
gas of the more massive system, which is inconsistent 
with our data. To reproduce the WL-SZ mass ratios we 
see in the data, the foreground cluster would have to 
be M 200 = 3 x IO^H^Mq (M 50 q = 2.2 x 10 14 /i" 1 M o ), 
which is different from the dynamical mass at 1.8<r; how- 
ever, taking the data together, this larger mass is no more 
or less favored than the central dynamical mass value, as 
the WL and dynamical masses and the convergence peak 
locations carry large uncertainties. 

In our second set of simulations, we model a plane- 
of-sky merger by placing at z = 0.28 two M200 = 
3 x 10 14 h- r M Q (M500 = 2.2 x 10 14 ft.^ 1 M Q ) clusters at 
the k peak locations, which roughly reproduces a single 
-W500 = 4.30 x 1O 14 /i _1 M0 cluster if naively summed. 
This configuration also reproduces the double peak in 
the convergence field and the suppression of the shear 
profiles within i?5oo- However, the mass estimates in 
this case agree well with the SZ derived mass, in con- 
trast to the foreground contaminant simulation and the 
real data. We would expect the SZ centroid to fall in 
between the two merging clusters of similar mass, which 
is consistent with the data. 

In summary, we cannot fully differentiate between the 
two scenarios of a foreground cluster or a merger, or 
quantify their exact effect on the WL and SZ mass es- 
timates, though both scenarios could explain the sup- 
pressed shear profile within i?5oo,sz- The spectroscopic 
data show strong evidence for a lower mass interloping 
cluster at z = 0.14, which could lead to large WL-SZ 
mass ratios as seen in the data; however, there is signif- 
icant uncertainty on the interloper's mass making this 
interpretation inconclusive. The merging scenario would 
be more consistent with the observed offset in the WL 
peaks from the SZ data and their coincidence with the 
BCGs at z — 0.28, but a merger alone would not explain 
the large WL-SZ mass ratio. Therefore, the data may be 
pointing to a combination of a plane-of-sky merger and a 
foreground interloper. In this section we have quantified 
the impact of either scenario on the WL and SZ mass 
estimates. Regardless, either scenario does not impact 
the goal of this work, which is to provide an unbiased 
test of the accuracy of the masses of an ensemble of SZ- 
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selected clusters. We therefore do not treat this cluster 
differently from the others in the mean ratio tests. 

7.3. Other Works Comparing WL and SZ-derived 
Masses 



AMI Consortium: Hurley- Walker et al. (2012)| mea- 



sured the mass of galaxy clusters with WL data from the 
CFHT and SZ data from the Arcminute Microkelvin Im- 
ager (AMI) . The four clusters for which both WL and SZ 
masses were successfully estimated exhibited a weighted- 
mean ratio of WL-to-SZ masses of 0.86 ±0.14, with total 
scatter of 14%. While this agrees with unity and with 
our baseline results at < la, such a direct comparison 
between their work and ours is not straightforward given 
the different overdensity used (A = 200), their use of a 
joint Bayesian analysis of the two data sets, and other 

d ifferences in analysis techniques. 

iPlanck Collaboration: Aghanim et al. (2012) showed 



that the amplitude of SZ signal derived from Planck 
data versus WL masses derived from Local Cluster Sub- 
structure Survey data obtained with the Subaru tele- 
scope disagrees at > la with results calibrated to X-ray 
data assuming hydrostatic equilibrium, but in the oppo- 
site sense than that predicted by non-thermal pressure 
support. That work proposes that systematic errors in 
the WL analysis could give rise to the discrepancy seen, 
and that further, careful study of systematics is needed 
to determine whether the difference is astrophysical in 
origin. The discrepancy is also se en in the scaling rela- 
tion result of Marrone et al. (2012 1, which used the same 
WL data but used SZ data from the SZ Array (SZA). 

8. CONCLUSION 

We have observed five galaxy clusters with the Mega- 
cam camera on the 6.5m Magellan Clay telescope, with 
the goal of measuring total masses with weak gravita- 
tional lensing and empirically testing th e accuracy of the 
SZ and X-r ay-based mass estimates of|Reichardt et al~ 



(2012 R12). Shear is estimated in deep r'-band data 



and additional g'- and i'-band data are used to calibrate 
photometry using the stellar color-color locus, as well as 
to remove cluster galaxies using color cuts. The source 
redshift distribution is estimated from CFHTLS-Deep 
photometric redshift catalogs (which do not overlap our 
fields) under the same photometric selections. 

We adopt three measures of total mass derivable from 
WL data: aperture masses, which are inherently two- 
dimensional quantities; spherical masses derived from 
NFW fits to shear data and evaluated at i?5oo,sz deter- 
mined from the SZ data; and the same NFW fits eval- 
uated at -Rsocwl as determined from the best-fit pro- 
files themselves. In all cases, only shear profiles from 
i?50o,sz to 12' are used. To make one-to-one compar- 
isons with the WL derived aperture masses, we com- 
pute the aperture mass equivalent statistic of Msoo.sz, 
which requires the assignment of an NFW profile using 
concentration-mass-redshift scaling relations from the 
literature. For the spherical mass comparison, a c—M—z 
relation is also adopted to fit models to the shear pro- 
files, and the resulting WL mass is compared directly to 
Af5oo,sz- Calibration tests are performed on mock cat- 
alogs based on ray-traced A^-body simulations. These 
show evidence for bias at levels consistent with a num- 
ber of previous works, under the assumption of perfect 



knowledge of cluster centers, galaxy redshifts, and galaxy 
shears. Under all methods, the mean ratio of WL to SZ 
masses is consistent with expectation of unity to within 
the statistical uncertainty: the mean aperture mass ra- 
tio is 1.04 ± 0.18, the mean spherical mass ratio using 
i?500,sz i s 1-07±0.18, and the mean spherical mass ratio 
using i?5oo,WL is 1.10 ±0.24. Consistency checks are per- 
formed, wherein we make more conservative selections on 
both the Clay-Megacam and CFHTLS-Deep photometric 
catalogs, and resulting mean ratios remain statistically 
consistent with unity and with the baseline results in all 
cases. 

Possible sources of systematic error are explored. The 
dominant sources are most likely the assumption about 
concentration and calibration to A'-body simulations. 
Different concentration scaling relations change mean 
aperture mass ratios by 2% to 6%, and a few percent 
more in the mean NFW mass ratios, depending on the 
relation adopted. Reducing this source of error will be 
key to unbiased and precise mass constraints, and this 
requires obtaining knowledge of the concentration of the 
cluster population that SPT detects. The bias levels with 
respect to Af-body simulations is consistent with levels 
seen in other works, at —6% to —13%. 

The assumed cluster center is another potential source 
of uncertainty. We show that using the BCG posi- 
tion give nearly exact agreement with baseline results 
in which the SZ position is used. While using the recon- 
structed convergence field peak gives the greatest devi- 
ations from the baseline result, this is probably because 
of noise in centering due to the large statistical uncer- 
tainties in determining convergence field peak locations. 
Other sources of systematic error, including shear bias, 
assumed cosmological parameters, and statistical uncer- 
tainty on the source redshift distribution, are subdomi- 
nant to these. 

For two of the five cluster systems in this work, we 
identify signs of possible merging activity and structures 
in the line-of-sight. We discuss these systems and simu- 
late the effects of the proposed contaminants. The sim- 
ulations show that both plane-of-sky merging activity 
and line-of-sight structures can induce multiple peaks in 
the convergence field and suppress shear profiles within 
-Rsoo.sz, but have different effects on the WL derived 
masses. The tests are not conclusive about the actual 
physical configuration of the clusters. The goal of this 
wor k is t o provide as unbiased a test of the accuracy of 
the |R12| mass estimates of SZ-selected clusters as possi- 
ble, so we do not treat these clusters differently in our 
analysis. 

In conclusion, we find statistical consistency between 
masses derived from WL data and those derived from 
SZ and X-ray data at the ~ 20% level. This represents 
first steps toward improved galaxy cluster mass estimates 
in the SPT survey. Improving the calibration of mass- 
observables is critical for exploiting the full statistical 
power of the SPT 2500 deg 2 survey data set in cosmolog- 
ical cluster abundance studies. 
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APPENDIX 
ADDITIONAL FIGURES 

The left panels of Figures [8j [10J pj [Mj and [16] show the SPT SZ detection significance maps. The SZ images 
subtend 20' on a side. The mapping between color and SZ significance £ spans the full range of SZ pixel values in 
the region of sky shown. The negative lobes are due to the filtering of the time-ordered data and the maps. Black 
contours correspond to the k reconstructed from the WL shear catalogs and denote values spaced by Ak = 1/30. The 
diamond symbols indicate the peak of the reconstructed n maps over the full field. Contours are dashed where values 
are negative, and solid where values are positive. The white X symbols indicate the centroid from the SZ significance 
map. 

The right panels of these figures show the Clay-Megacam images. A false-color composite is presented, with irg 
mapped to the RGB channels, respectively. The same k contours from the left panels are shown in cyan. White contours 
correspond to the SZ significance data from the left panels, spaced at significance values of (—8, —4, —2, 0, 2, 4, 8). The 
horizontal white lines enclose the BCG of each cluster. 
Shear and aperture mass profiles are shown in Figures [9} 1 1 1 1 [13} 1 1 5 1 and|17J The shaded region denotes radii in terior 



to i?5oo,sz- The left panels show the binned shear data wlTh DesWit NFW profiles assuming concentrations of D08 



These models are fit at radii 9\ < 9 < 9 2 , where 9\ corresponds to i?soo,SZ- The right panels show the aperture mass 
profiles within circular apertures of radius 9\. For all data, 9 2 has been fixed to 12'. The cyan shaded region denotes 
the ±lcr s tat excursion in the aperture mass at a given 9\. Cluster masses inferred from SZ data are shown a s data 
points with error bars. These masses have been projected assuming spherically symmetric, three dimensional NFW 
profiles with |D08| concentrations, and then filtered using the same kernel as the WL aperture mass statistic at the 
given 9\. 



20 



F. W. High et al. 





Fig. 8. — SZ, optical, and k data for SPT-CL J0516-5430. See the Appendix for a description. 
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Fig. 9. — Shear and aperture mass profiles of SPT-CL J0516-5430. See the Appendix for a description. 
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Fig. 11. — Shear and aperture mass profiles of SPT-CL J2022-6323. See the Appendix for a description. 
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Fig. 12. — SZ, optical, and n data for SPT-CL J2030-5638. See the Appendix for a description. 
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Fig. 13. — Shear and aperture mass profiles of SPT-CL J2030-5638. See the Appendix for a description. 
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Fig. 14. — SZ, optical, and ft data for SPT-CL J2032-5627. See the Appendix for a description. 
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Fig. 15. — Shear and aperture mass profiles of SPT-CL J2032-5627. See the Appendix for a description. 
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Fig. 16. — SZ, optical, and k data for SPT-CL J2135-5726. See the Appendix for a description. 




Fig. 17. — Shear and aperture mass profiles of SPT-CL J2135-5726. See the Appendix for a description. 



